Phase diagram of turbulence in superfluid 3 He-B 

A.P. Finne*, S. Boldarev*t, V.B. Eltsov* t, and M. Krusius* 

* Low Temperature Laboratory, Helsinki University of Technology, 

P.O. Box 2200, FIN-02015 HUT, Finland 
< Kapitza Institute for Physical Problems, 119334 Moscow, Russia 



In superfluid 3 He-B mutual-friction damping of vortex-line motion decreases 
roughly exponentially with temperature. We record as a function of temper- 
ature and pressure the transition from regular vortex motion at high temper- 
atures to turbulence at low temperatures. The measurements are performed 
with non-invasive NMR techniques, by injecting vortex loops into a long 
column in vortex-free rotation. The results display the phase diagram of 
turbulence at high flow velocities where the transition from regular to tur- 
bulent dynamics is velocity independent. At the three measured pressures 
1 0. 2, 29. 0, and 34 bar, the transition is centered at 0. 52 - 0. 59 T c and has 
a narrow width of 0.06 T c while at zero pressure turbulence is not observed 
above 0.45 T c . 

PACS numbers: 47.37, 67.40, 67.57 
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1. INTRODUCTION 

According to the original Landau picture of superfluidity superflow is ir- 
rotational: V x v s = 0. This condition severely restricts the possible motions 
of the superfluid: When placed in a rotating container, it cannot participate 
in rotation. Later Onsager and Feynman suggested the modern view of su- 
perfluid dynamics: By the creation of quantized vortex lines the superfluid is 
able to mimic arbitrary complex flows on length scales which are large com- 
pared to the inter-vortex spacing. It was discovered more than 40 years ago 
that if superfluid 4 He is driven sufficiently fast then a complex vortex tan- 
gle appears. This tangle produces apparently chaotic time-dependent flow 
at different length scales. Such motion was called superfluid or quantum 
turbulence.^ 

Studies of turbulence in a different superfluid, the B phase of super- 
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fluid 3 He, have started only recently!^ From the hydrodynamics point of 
view 3 He-B differs from 4 He in two important respects. First, the viscosity 
of the normal component in 3 He-B is large compared to both the (sample 
size) x velocity in a typical experiment and to the circulation quantum. This 
means that the normal component of 3 He-B is always in well-defined ex- 
ternally imposed laminar motion, which considerably simplifies the analysis 
of the experimental results. Second, the damping of the superfluid motion 
varies significantly in 3 He-B as a function of temperature. Such damping 
comes from the interaction of the thermal quasiparticles with the cores of 
the quantized vortices (so-called mutual friction). The magnitude of the 
frictiorPI changes from large values (typical for superconductors) in the limit 
T — > T c to small values (typical for 4 He-II) at temperatures T < 0.5 T c . 

In a recent experiment we discovered that there is a sharp transition 
in the character of superfluid motion in 3 He-B as a function of tempera- 
ture at P = 29 bar pressure. When the laminar vortex-free superflow is 
prepared by rotation of the cylindrical sample and a few seed vortex loops 
are injected into it, then these loops expand in a regular manner, with their 
number conserved at T > 0.6 T c or, proliferate to a turbulent vortex tangle 
at T < 0.6 T c . (Eventually the turbulent tangle polarizes to mimic the global 
rotating flow of the normal component which compensates the driving force 
and the turbulence decays to a cluster of rectilinear lines containing ~ 10 3 
vortices.) 

The transition to turbulence is controlled by the mutual friction damp- 
ing. It can be understood by analyzing the dynamic equation for the coarse- 
grained superfluid velocity v s which is averaged over volumes containing 
many vortex lines P 

dv ( v \ 

+ V( /x+ y J = v s x u + a'u x (v s -v n ) + au x [u x (v s -v n )]. (1) 

Here v Q is the flow imposed on the normal component, u = V x v s , u> is a 
unit vector in the direction of u, \i is the chemical potential, a and a' are 
the dissipative and reactive mutual friction coefficients, respectively. The 
vortex line tension is neglected. Applying simple dimensional comparison of 
the magnitudes of the dissipative term (containing a) to the inertial terms in 
Eq. I|T|). similar to how one introduces the Reynolds number in the classical 
Navier-Stockes equation, one arrives to the conclusion*^ that the character 
of the flow is controlled by the ratio of the mutual friction coefficients q = 
a/(l — a'): Turbulence occurs only if q < q c ~ 1. The transition from 
laminar to turbulent dynamics becomes velocity-independent in the range 
of validity of Eq. Q (which includes the requirement for a sufficiently large 
flow velocity). This is seen from the scale invariance of Eq. If v s (r,t) 
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is the solution of this equation (and of the continuity equation div v s = 0) 
for a given imposed flow v n , then Xv s (r, Xt) is the solution for the scaled 
imposed flow Xv n . 

In this report we extend the investigation of the transition between 
regular and turbulent vortex dynamics in 3 He-B to the pressure range from 
zero up to the solidification pressure and examine the question whether the 
transition to turbulence is controlled by a single universal value q c at all 
conditions. 

2. EXPERIMENT 

A detailed account of the experimental techniques, including the injec- 
tion of the seed vortex loops and the NMR detection, is given in Ref . EJ Here 
we present a brief summary. The 3 He sample is contained in a cylindrical 
container with radius R = 3 mm and height 11 cm. At mid- height there is a 
barrier magnet which provides a magnetic field to stabilize the A phase at 
all temperatures. The lower and upper sections of the sample remain in the 
B phase. Two NMR detection coils are installed in these sections close to 
the top and bottom ends of the sample. From the NMR measurements the 
number of vortex lines within the coil can be determined. 

In the absence of the A phase the B phase in our experiments remains 
vortex-free in rotation up to some container-specific critical angular velocity. 
Its value is not in good control but for the containers used in the present 
work it exceeds 2rad/s at T < 0.8 T c . When the A phase is present in the 
middle of the sample the seed vortex loops are injected into the B phase at 
a well defined and reproducible as a result of the Kelvin-Helmholtz (KH) 
instability of the AB interface. 7 This velocity is a well-understood smooth 
function of T, P and of the current in the barrier magnet lb- The circu- 
lation carried by the initial B-phase vortices comes from the vortex layer 
which covers the AB interface. In the layer vorticity is arranged in a regu- 
lar manner, with vortices aligned radially and parallel to the AB interfaceP 
Thus one may expect that the configuration of the seed vortices is reasonably 
well reproducible with one end of each vortex sticking out of the AB inter- 
face and the other end perpendicular to the outer wall of the sample. Our 
measurements indicate that the size of a seed vortex is about 0.4 mmP The 
injection parameter, which is not in good control, is the number of the seed 
vortices, produced by one instability event. This number has a rather wide 
distribution in the range 3-30 vorticesj^ with the average value of about 10. 
The reason for such a wide variation is presently not known. 

The result from the expansion of the seed vortices along the B-phase 
column is observed when they reach the NMR coils. A few vortices (in the 
case of regular expansion) and an array with many vortex lines (in the case 
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Fig. 1. Transition between regular (open squares) and turbulent (filled 
squares) vortex dynamics as a function of temperature at three pressures. 
The dashed vertical lines represent the average values T t of the distributions 
in Fig. H 



of turbulent expansion) are easily distinguished at T > 0.4 T c . The two B- 
phase sections of the sample present two independent experimental volumes. 
Their simultaneous measurement improves the statistics of the results. 



3. RESULTS 

Figs. HE] summarize the transition from regular to turbulent dynamics 
with decreasing temperature and trace the transition as a function of Q, 
T, and P. Here is the critical velocity for the KH instability, at which 
vortices have been injected in the flow. In most cases the B phase has been 
completely vortex-free before the injection and thus the flow in the sample is 
v n — v s = fi x r. The data have been measured along continuous trajectories, 
either by scanning temperature at constant barrier current lb or by scanning 
lb at constant T. These trajectories are not emphasized in the plots, instead 
the data have been classified in regular (open symbols) or turbulent (filled 
symbols) events, to highlight the transition separating them. 

On the basis of the data in Fig. ^ (which is plotted in a different way 
in the later figures) we conclude that the transition is independent of the 
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Fig. 2. Probability for a vortex expansion event to be regular (open squares) 
or turbulent (filled squares) as a function of temperature at three pressures. 
The probabilities are calculated by arranging events from Fig. ^ as a func- 
tion of temperature in bins of 0.01 T c width at the three pressures. The 
continuous curves are fitted normal cumulative distribution functions with 
the average transition temperature T t and dispersion at, as given in each 
panel. 

velocity SI, if compared with the width of the transition. This agrees with 
expectations based on the scale invariance of the equation Q . This conclu- 
sion is valid for high velocities, while low velocities much below the present 
range of ~ lrad/s should be checked separately. 

Other features demonstrated by Fig. ^ are that the transition has a 
certain width and that the transition temperature might depend on pressure. 
To quantify these features the probabilities to oberve regular and turbulent 
events are plotted in Fig.|2]as a function of temperature for the 3 pressures. 
Each of these distributions is fitted with a normal cumulative distribution 
function with average value Tt and dispersion at as fitting parameters. This 
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function has been chosen as a simple way to obtain numerical values for the 
average transition temperature and its width. We do not claim that the 
usage of the Gaussian distribution has a solid physical reason. 

For all 3 pressures the transition half-width is found to be around 
<Tt ~ 0.03 T c . We believe that fluctuations in the injection process are the 
main cause for the widths in Fig. It appears probable that at higher tem- 
peratures only some of the initial configurations of vortices (say, the ones 
with a larger number of initial loops) will evolve to a turbulent network. As 
the temperature decreases more of the different initial configurations become 
unstable towards turbulence. This picture is supported by comparison with 
experiments on turbulence initiated by neutron absorption in 3 He-BP^KH 
injection which we use in this work produces almost always more than 3 vor- 
tex loops, while in the neutron absorption the production of 1-2 vortex loops 
is the most common case. As has been observed in Ref. I10( the transition 
to turbulence becomes wider towards lower temperatures if it is initiated 
by neutron absorption: Even at 0.45 T c regular events are still rather com- 
mon. Numerical simulations of vortex dynamic*^ also demonstrate that 
at intermediate temperatures a single vortex loop in the rotating container 
may initiate turbulence or may not, depending on the initial configuration. 
Also in the flow of classical liquids in a pipe the transition from laminar to 
turbulent flow depends on the initial perturbation: The critical Reynolds 
number has been found to scale inversely proportional to the perturbation 
amplitude.^ 

The dependence of the average transition temperature T t on pressure is 
shown in Fig. El The data from Fig. |2] have been augmented by the measure- 
ments at P = 0. At zero pressure we have not observed the turbulence at 
T > 0.45 T c at rotation velocities 0.5 - 0.7rad/s. On the other hand, turbu- 
lent vortex tangles have been observed at T < 0.2 T c using vibrating wiresP 
Combining all the data, we conclude that there is a pressure dependence of 
the transition to turbulence on the T/T c scale. 

When inspecting Fig. |3] one should keep in mind that from the three 
sets of measurements in Fig. ^ more effort was invested in the tempera- 
ture calibration and data taking of the runs at 29.0 and 10.2 bar pressures, 
while the 34 bar measurements were more qualitative in character. Thus, 
a non- monotonous dependence of T t /T c on pressure cannot be claimed. It 
may simply be an artefact of the small amount of data measured at 34 bar. 
However, we note here an interesting phenomenon, which might affect the 
transition to turbulence at high pressures. As indicated by the dashed line in 
Fig. |31 there is a first-order transition^ in the vortex core structure between 
a non-axisymmetric core at low T, P and an axisymmetric core at high T, 
P. The structure of the core affects the mutual friction parameters and the 
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Fig. 3. Phase diagram of turbulence in 3 He-B. The points for 10, 29 and 34 
bar represent the middle of the transition T t from regular dynamics at higher 
temperatures to turbulent dynamics at lower temperatures, as determined 
from the distributions in Fig. [21 The error bars indicate the width of the 
transition o" t as determined from these distributions (and do not indicate 
the uncertainty in temperature). At zero pressure the error bar indicates 
the unexplored region where the transition to turbulent dynamics should 
occur. The solid line is a guide for the eye. The dashed line shows the 
transition in vortex core structure^ which leads to a discontinuity in the 
mutual friction coefficients.^ 

dissipation associated with a non-axisymmetric vortex is larger. 4 It is there- 
fore possible that there is a kink in the T t (P) dependence at the transition 
between the two core structures. 

What is the reason for the pressure dependence of the transition to 
turbulence? Eq. Q suggests that the transition is controlled by the ratio 
of the mutual friction parameters, q = a/(l — a'). Thus it is instructive to 
plot our data as a function of q. Unfortunately, the precise measurements 
of q(T, P) are available only for two pressures, 10 and 29 barP Our data 
for these pressures from Fig. ^ have been repeated as a function of q in 
Fig. IU Generally speaking, the transition takes place when q ~ 1. However, 
a critical value of q is not universal: the transition moves to larger q as 
a function of pressure. Can we be sure that the difference in q values for 
the transitions in Fig. ^ is experimentally significant? First we note that the 
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Fig. 4. Transition between regular (open squares) and turbulent (filled 
squares) vortex dynamics as a function of the ratio of the mutual friction 
coefficients q = a/(l — a') at the two pressures where q(T, P) was measured 
in Ref . H 



mutual friction coefficient a can be measured in our experiment in situ using 
the time-of-flight technique described in Refs.EHandlEl Such measurements 
give values of a which agree with Ref.|l]both at 29 and 10 bar pressure. Thus 
there is no systematic error, for example, in the temperature scales. The 
scatter of the data points in Ref. 4 gives an uncertainty in q of about ±0.04 at 
P = lObar. For P = 29 bar the situation is more difficult to estimate, since 
there exists an additional scatter due to the presence of two different types 
of vortices around 0.6 T c . The possible range of q values at the transition 
temperature is from q ~ 1.13 for the high-temperature axisymmetric vortex 
to q ~ 1.72 for the low-temperature non-axisymmetric vortex (the same type 
which occurs at 10 bar). Thus, even with these uncertainties, the critical 
values of q for the two pressures in Fig. 0] remain well separated. 

Thus, whether the vortex dynamics is turbulent or regular in charac- 
ter, indeed strongly depends on the magnitude of mutual friction damping. 
However the simple classification on the basis of the value of q, as suggested 
in Ref. 121 agrees with the experiment only as an order of magnitude esti- 
mation. To explain the pressure dependence of the transition to turbulence 
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presented here more experimental and theoretical work is required. At the 
moment we can only speculate about possible reasons. 

First, the analysis of Eq. may not be completely correct. This equa- 
tion includes two parameters, a and a' , and the transition may depend on 
their values in a more complicated way than the simple ratio q. Second, the 
equation itself may not be fully applicable to the experimental situation. In 
particular, in the process of coarse- graining the mutual friction force over 
volumes containing many vortex lines the values of the mutual friction coeffi- 
cients may become renormalized from their single- vortex values measured in 
Ref.QJ Next, the properties of the injection of the seed loops may depend on 
pressure in such a way that the initial configurations of vortices become less 
favourable for initiating the turbulence with decreasing pressure. However, 
in the latter case one would expect that the width of the transition will also 
increase, similar to the turbulence induced by neutron absorption,^ which 
is not observed in Fig. |2J 

At present it is assumed that one vital precondition is required to start 
turbulence: it is the Kelvin-wave instability^ of a vortex line. It provides 
the mechanism by which quantized vorticity starts to multiply even from one 
single vortex loop which is injected in the rotating vortex-free counterflow. 
The Kelvin-wave instability switches on when some section of the vortex loop 
becomes oriented along the flow. In this part of the vortex filament a helical 
Kelvin wave starts to grow in amplitude. The wave develops into loops which 
reconnect to form separated vortex rings. These rings provide the beginning 
of the evolving network. The triggering and development of the Kelvin- wave 
instability is significantly affected by the vortex line tension, which has so 
far been neglected here. The tension depends, although only logarithmically, 
on the vortex core size, which decreases in 3 He-B by almost an order of 
magnitude from zero to melting pressure. This might also influence the 
transition to turbulence. 

The final complication, which we note here, concerns the role of the 
magnetic field. Close to the AB interface, where the seed vortex loops are 
injected, the magnetic field is high and distorts the superfluid gap in the B 
phase as well as the structure of the vortex core. These are the parameters 
which directly determine the magnitude of the mutual friction coefficients.^ 
If the turbulence is initiated in the immediate vicinity of the original position 
of the seed loops then the transition to turbulence might also be affected by 
the pressure and temperature dependent magnetic field Hp&. The measure- 
ments with neutrons^ are free from these complications but the transition 
in this case is wide, as has been explained above, and the presently available 
results do not allow to pinpoint the critical q value with enough accuracy to 
judge its pressure dependence. 
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4. CONCLUSIONS 

We have described the first measurements on the phase diagram of tur- 
bulence in 3 He-B: (i) A wide regime of mutual friction damping has been 
examined as a function of temperature, which is not easily accessible in the 
case of 4 He-II. Mutual friction proves to control the character of vortex mo- 
tion in superfluid hydrodynamics: when a>/(l — af) > 1 the number of vortex 
lines is a conserved quantity, while in the opposite case, a/(l — a') < 1, the 
vortex number rapidly proliferates through turbulence. It is known from 
4 He-II that mutual friction can help to sustain turbulence by blowing up 
Kelvin-wave perturbations on vortex loops and it also can damp turbulent 
motion at certain length scalesP^^From 3 He-B we now find that high mutual 
friction can totally suppress turbulence, (ii) The transition to turbulence is 
not marked by a unique value of the mutual friction ratio q = a/(l — a'); 
with increasing pressure the transition moves to higher q. (iii) Turbulence is 
prominently concentrated in the low-temperature, high-velocity, and high- 
pressure corner of the 3 He-B phase diagram. At higher flow velocities the 
transition between regular and turbulent dynamics is velocity independent, 
the width of the transition regime is narrow and centered in the temperature 
interval 0.5 - 0.6 T c , depending on pressure. 
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